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Abstract 

We study the inviscid limit problem of the incompressible flows in the presence of 
both impermeable regular boundaries and a hypersurface transversal to the boundary 
across which the inviscid flow has a discontinuity jump. In the former case, boundary 
layers have been introduced by Prandtl as correctors near the boundary between the 
inviscid and viscous flows. In the latter case, the viscosity smoothes out the discontinuity 
jump by creating a transition layer which has the same amplitude and thickness as the 
Prandtl layer. In the neighborhood of the intersection of the impermeable boundary and 
of the hypersurface, interactions between the boundary and the transition layers must 
then be considered. In this paper, we initiate a mathematical study of this interaction 
and carry out a strong convergence in the inviscid limit for the case of the plane parallel 
flows introduced by Di Perna and Majda in [5]. 

1 Introduction 

In this paper we are interested in the behavior of the incompressible Navier-Stokes flow 
when the viscosity is small. This so-called inviscid limit problem is particularly difficult 
when the flows is contained in a domain limited by impermeable walls. In the standard case 
of a half-space, the problem reads as follows: 



dtv^ + (v' • V)v' + Vp' = eAv' 
div v' = 0. 



(1.1) 



Here, x = {x,y,z) is in M x M x (0, +oo), the velocity v*^ = {u'' , v'' , w'') is in M^, p*^ the 
pressure and e > is the viscosity parameter. 

The equation (II. ip is imposed with the classical no-slip boundary condition: 

vf =0. (1.2) 

|z=0 ^ ^ 
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Considering the problem (jl.ip - (|1.2p in the hmit e — )• 0, one may hope to recover the Euler 
flow: the equation (jl.ip with e = 0, for which the natural condition on the boundary {z = 0} 
is 

"'to = 0- (1-3) 

Due to the difi^erence (or rather, loss) of boundary conditions, it is common in the 
limit to add a boundary corrector or the so-called Prandtl layer. This formal procedure 
was introduced by Prandtl in 1904, and it remains a challenging mathematical problem to 
circumvent the validity of this theory. 

Yet, some positive answers have been given in the setting of analytic flows in two di- 
mensions by Caflish and Sammartino in [S] and improved by Cannonne, Lombardo and 
Sammartino in [lOj for inviscid flows that are analytic with respect to the tangential vari- 
ables. On the other hand, when the smoothness with respect to the tangential variables is 
limited, the Prandtl layer have been shown to be unstable; see for example the papers by 
Grenier [4], Gerard- Varet and Dormy [3j, Guo and Nguyen [5]. Note that these papers also 
concern the 2d case. 

Here, we propose to study the inviscid limit problem of a viscous incompressible Navier- 
Stokes flow in presence of both a solid boundary and of a transversal discontinuity hypersur- 
face in the limiting inviscid flow. The full problem is currently out of reach. In particular, 
jump discontinuity across a hypersurface is also a rather unstable pattern for the incompress- 
ible Euler equations, because of the Kelvin-Helmhotz instabilities. Nevertheless, when the 
inviscid theory is successful to provide some Euler solutions with some jump discontinuity 
across a hypersurface, it is expected that the extra viscosity in the Navier-Stokes solutions 
smoothes out the discontinuity into a transition layer which can be basically thought as 
a transmission version of the Prandtl layers. Here again, positive results are known in an 
analytic framework, in 2d, see [8]. 

However, since this hypersurface is assumed here to be transverse to the boundary, one 
cannot relies on the previous results based upon the analyticity in the transversal variables 
to study the interactions between the boundary layer and the transition layer. We will 
therefore study the layers interactions in Sobolev spaces. 

In this paper, we will restrict our study to a simple setting of three-dimensional incom- 
pressible flows: the plane-parallel flows. They were introduced by DiPerna and Majda in 
[2j in order to prove that the Euler equations are not closed under weak limits (in three 
spatial dimensions). These flows have also been used as basic flows for the Euler equations 
by Yudovich in to investigate stability issues and recently by Bardos and Titi in [Ij to 
investigate several longstanding questions including the minimal regularity needed to have 
well-posedness results, localization of vortex sheets on surfaces, and the energy conservation 
for the Euler equations. 
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Precisely, a plane-parallel solution is of the form: 




^^'{t,x,y,z)= ^'{t,x,z)\ . (1.4) 
Then, the Navier-Stokes system (jl.ip depletes into 

= ^^^^.^ 

with ]f = 0. It is thus a pressureless flow. Observe that a vector field of the form ()1.4p is 
divergence free. On the other hand the boundary conditions p.2p now read 

iu\v%^^ = 0, (1.6) 

as the Dirichlet condition for the third component is automatically satisfied for flows of the 
form ()1.4p . The system ()1.5p - ()1.6p is now quite simple: the first equation in ()1.5p is a one 
dimensional heat equation whereas the second one is a two dimensional transport-diffusion 
equation, and for both we prescribe homogeneous Dirichlet conditions. 

On the other hand, the Euler system, the equations (II. Ij) with e = 0, depletes into 

dtu^ = 

dtv^ + u^d^v^ = 0. ^^'^^ 
Therefore the solution starting from the initial data 

vo(x,z) = I vo{x,z) 







is simply given by the formula 



uoiz) 

v°(t,x,z)= \ vo{x - tuo{z),z) I . (1.^ 

This holds true in a quite general setting, but let us be formal for a few more lines. For 
instance let us think that the function vq is smooth for a while, so that there is no doubt to 
have about the meaning of the formula (jl.Sp nor about the fact that it solves the depleted 
Euler equations (II. 7p . We want to focus here first on the issue of the boundary conditions. 
In particular, note that no boundary conditions are needed to be prescribed for the system 
(jl.7p . since any solution of the form ()1.4p already satisfies the condition ()1.3p . On the other 
hand, if the initial data vq does not vanish on the boundary z = then neither does the 
corresponding solution v'^ given by (jl.Sp for positive times. As a consequence, v'' does not 
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satisfy the condition (jl.6p and therefore cannot be a good approximation, say in L°°, of 
any smooth solution v*^ of the system (|1.5p - (jl.6p . Yet Prandtl's theory predicts that the 
system (|1.5p - (|1.6p admit some solutions v*^ which have the following asymptotic expansion 
as e — )■ 0: 

v^(t, X, z) ~ vO(t, X, z) + wp{t, X, (1.9) 

Above the profile vp{t,x, Z) describes a Prandtl boundary layer correction. In particular 
it satisfies vp(t, x,Z) — t- when Z — )• +oo, so that this term really matters only in a layer 
of thickness i/e near the boundary {z = 0}, and also satisfies v''(t, x, 0) + vp(t, x, 0) = 0, so 
that the functions in the right hand side of (|1.9p satisfies the boundary conditions (II. 6p . The 
validity of this asymptotic expansion has been verified in a recent paper of Mazzucato, Niu 
and Wang [7] for regular initial data vq. In particular it follows easily from their analysis 
that for any regular initial data vq, there exists a sequence of smooth solutions v*^ of the 
system (|1.5p - (|1.6p . with some initial data conveniently chosen, such that v*^ converges to 
strongly in the topology. 

Here, as mentioned previously, we are interested in the case where vq has a jump of 
discontinuity across a hypersurface. More precisely we assume that uq is smooth and that 
Vq is piecewise smooth with a jump of discontinuity across the hypersurface {x = 0}: 

N]l^o •= ^™ voix,z)- lim vo{x,z)^0. (1.10) 

x— >0+ a;— >0~ 

We assume for simplicity that there is no jump of the normal derivative of vq across the 
hypersurface {x = 0}, that is 

[dxVo]\^=,=0. (1.11) 

Then it can be easily seen on the formula (jl.Sp that the corresponding Euler solution 
is piecewise smooth with a jump of discontinuity across the hypersurface given by the 
equation {"^^{t,x, z) = 0}, with 

^^{t, x, z) := X — 'ilj{t, z), 'ip{t, z) := tuo{z), 

Moreover taking the derivative with respect to x of the both sides of Formula (jl.Sp yields 
that there is no jump of the normal derivative of vq across the hypersurface {"^^{t, x, z) =0}. 

Such a pattern cannot hold anymore for any reasonable solutions v*^ of the depleted 
Navier-Stokes equations (jl.5p - (11.6p : the viscosity smoothes out this jump of discontinuity 
into a transition layer near the hypersurface = 0}. In particular v*^ and its normal 
derivative must be continuous across the hypersurface = 0}: 

[v1|,o.o=0a^d[a.v^]|^,^^=0. (1.12) 

Following Prandtl's ideas it is natural to introduce a corrector 

^°(t,x,z) 
^rKH{t, ^ ,z) 
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where the profile VKH{t,x,X) satisfied 

[Vkh]\x=0 = -Nx=0 /-I -ION 

[dxVKH]u^, = 0, ^'■'-'^ 

and Vkh — s- as X — t- ±00. This strategy can be seen as a transmission counterpart of the 
introduction of the boundary layer vp previously mentionned. Actually, if the fluid domain 
was not limited by the boundary {z = 0} one could then adapt the analysis of [7] to justify 
the existence of some solutions v*^ of (ll.5p - p.6p which admits an expansion of the form 

V [t, X, Z) ~ V^'^t, X, Z) + WKH{t, j= ,Z). 

Yet there is no reason for which the transition layer ^kh should satisfy the boundary 
condition at 2; = 0, nor for which the boundary layer vp should take care of the jump 
condition across {"^^ = 0}. It is precisely our point to understand how to deal with both 
layers. 

Our result is the following. 

Theorem 1.1. Let 1 < p <2 and let 

Uo{z) E i72(0,+oo), 
VQ^+{x, z) E VF2'P([0, +00) X (0, +00)), (1.14) 
vo,.{x,z) E l^2,P((_oo,0] X (0,+oo)), 

and 

t;oELP(Mx(0,+oo)), w^th ^o(x, z) := | ^ > [j; (1.15) 

Assume that vq satisfies the jump conditions (jl.lOp and (jl.lip . Let us consider v*^ given 
by the formula (11. Sp . which for any T > is a distributional solution of the depleted Euler 
equations (|1.7p with V||_^ = vq. 

Then, there exist some smooth solutions v*^ := {u''{t, z),v^{t, x, z)) of the depleted Navier- 
Stokes equations (ll.5p - (ll.6p such that as e — t- 0, there holds the convergence 

,0 



in L°°{0,T;L\R+) X LP (RxR+)). (1.16) 



Here VF^'^ denotes the usual Sobolev space of order 2 associated to the Lebesgue space 
LP and Ll"^ denotes the special case Ll"^ := W"^'"^. 

Let us end our Introduction by giving here a few comments. 
^Here the subscript KH holds for Kelvin-Helmhotz 
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First, observe that in the statement of Theorem 11.11 the initial data of is not pre- 
scribed. In the proof, we will explicitly choose them in a convenient way; in particular, it 
allows the boundary and transmission layers to be initially specified. This could perhaps 
seem a little bit unusual at first glance, but it is in fact only technical for our convenient 
formulation of the main result. However this way to formulate our results avoids some extra 
considerations regarding forcing terms and/or initial layers which do not seem essential for 
our purpose in the present paper. 

Finally, let us mention that we are unable to include the case p = 2 or any p > 2 in 
Theorem II. 1[ We will explain why in Remark 13.71 



2 Straightened interface 

To fix the interface, we introduce the following change of variable: 

X := X — '4'{t, z), where '4'{t, z) := tuQ{z). 

In these coordinates, the discontinuity interface is given by the equation x = 0. In what 
follows, we drop the tilde in x. The system (II. Sp now reads 

(2.1) 

dtv' + {u' - uo)d,v' = eAtv', 

with 

At := dl + {d, - d^^d^f = (1 + m\'')dl - 2d,^dl, - dl^d^ + dl 

The boundary conditions (jl.6p do not change: 

(u^^;^)|^^„ = 0. (2.2) 



We are looking for some functions v!" and v'^ which satisfy the equations (j2.ip on both 
quadrants (x, z) G (0, +oo) x (0, -|-cxd) and (x^z) G (— oo,0) x (0, +00) with the interface 
conditions: 

[v^]|_„ = Oand[a,.v^]|_„=0, (2.3) 

which correspond to the conditions (I1.12P in the new variables. Now, since v!" does not 
depend on x, the conditions ()2.3p reduce to 

b1u.o = Oand[a,t;^]|_„ = 0. (2.4) 

Note that if v!" and are distributional solutions of (12. ip on both quadrants (x, z) G 
(0,-|-oo) X (0,-|-oo) and (x^z) G (— oo,0) x (0, +00) and satisfy the previous interface con- 
ditions then they are distributional solutions of (12. ip on the whole half-space M x (0, -|-oo). 

In the limit case e = 0, the situation is now particularly simple: in the new coordinates 
the solution v*^ is stationary 

' ^to(^) 

v"(t,x,z) = I t>o(x,z) I . (2.5) 
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Now, to prove Theorem II .11 it suffices to prove that there exist some functions w*^ and v*^ 
which satisfy the equations (j2.ip on both quadrants, satisfy the conditions (I2.2p and (|2.4p 
and converge, as e 0, to v° given by ([23]) in L°°(0,T; L2(M+) x LP{R x 



3 Asymptotic expansions 

Let us now describe our strategy. We are going to construct a family of functions of the 
form 

Z X X Z \ J 

Vlpp{t, X, Z) = Vq{x, z) + Vp{t, X, + VKH{t, Z) + ^^(t, 

Ve 

which satisfy approximatively (j2.ip on both quadrants (in a sense that we will precise in 
the sequel), and which satisfy the conditions ()2.2p and (]2.4p . 

In (|3.ip . (uc'Wo) are the functions given by ()2.5p . The other functions will be defined in 
the sequel. For instance, (?7p, Vp) will be a depleted Prandtl layer near the boundary, Vku 
will be a transmission layer near the discontinuity interface, and 14 will aim at describing 
the behavior of the boundary layers interaction. In what follows, we will use, as in the 
introduction, the capitalized variables X^Z to refer to xj ^^zj correspondingly. 

Then we will prove that there exists a family of functions (v!^ ^v^) close to (u^pp,t>^pp) 
which exactly satisfy (|2.ip on both quadrants and G (— oo,0] x (0,+oo), and which 

still satisfy the conditions ()2.2p and (|2.4p . 

It will remain to prove that this family (ti*^, u*^) converges to (no, vq) in L°^(0, T; L^(M4_) x 
LP(M X M_(.)), for 1 < p < 2, to conclude the proof of Theorem ll.il 



3.1 Construction of the approximated solution 
3.1.1 Construction of Up 

We start by defining the function C/p, which aims at compensating the non- vanishing value 
of at z = 0. On the other hand, we want this correction to be localized near the boundary 
z = 0. We will therefore require Up to satisfy 

f7p(t,0) = -tto(O), lim C/p(t,Z)=0. (3.2) 

Now if we put the Ansatz (13. ip into the system (12. ip and match the order in e, we get from 
the equation for v!" the following equation for the profile Up{t, Z): 

dtUp = dlUp. (3.3) 

We choose for Up the initial value 

Up\^^,{Z) = -uoifi)e-^. (3.4) 
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By Duhamel's principle, the solution Up{t,Z) of (|3.2p - (|3.3p - (|3.4p satisfies 

Up{t,Z) =-uo{0)e'^ -uo{0) / g{t- s,Z;Z')e-^'dZ'ds 

Jo Jo 

where Q{t, Z; Z') denotes the one-dimensional heat kernel on the half-line: 



g{t, Z; Z') := G{t, Z-Z')- G{t, Z + Z') 



G{t,Z) :-- 



'A-Kt 



-.e 4t . 



(3.5) 



Now by using the standard convolution inequality: < we easily deduce 

that, for any p > 1, 



\Up\\l°°{0,T;Lp{^+)) < Cq\uq{Q)\ 



1 + 



\\G{t,.)\\LA\e 



LP 



ds 



< Co\uoiO)\, 



for some positive constant Cq that depends on p and T. Here, we used the fact that 
= 1- Similarly, using the fact that ||5^G(t, ~ we obtain 



^zf^p||L°°(0,T: 



;LP(M+)) < C'o|no(0)| 1-F sup / (t-s) ^^^ds 

L 0<t<TJo 



which is again bounded by Cq \ uq (0) | . 

That is, we obtain the following lemma: 

Lemma 3.1. There exists a unique solution Up to the problem (|3.2p - (|3.3p - (j3.4p on [0,T] x 
R+, for any T > 0. Furthermore, for any p > 1, there is some positive constant Cq that 
depends on p and T such that 



\\Up\\l°°{o,T;W^-p{r+)) < Co\uo( 



(3.6) 



3.1.2 Construction of Vp 

For Vp, the situation is the same as that for Up, other than the fact that Vp also depends 
on the variable x. However, x only appears as a harmless parameter. More precisely, by 
plugging the Ansatz (j3.ip into the system (|2.ip and match the order in e, we then get from 
the equation for the profile equation for Vp{t,x, Z): 



dtVp = dlVp, Vp\^^^ = -7;o(x,z)|^^^, ^lim^^P = 0. (3.7) 

Once again, we choose an initial data compatible with the boundary condition, for instance 

Vp\^^, = -vo{x,Q)e-^. (3.8) 



Then as was the case for Up, there exists a unique solution Vp of (13.7p - (|3.8p satisfying the 
Duhamel principle: 



ft f + OO 

/ / Git- s, Z; Z')e-^' dZ'ds, 
Jo Jo 



Vp{t,x,Z) = - vo{x,Qi)e ^-uo(x,0) 
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where Q{t — s,Z;Z') is the heat kernel defined as in (j3.5p . It is clear from this integral 
representation for Vp that the only dependence on x is due to vq(x,0). Thus, we easily 
obtain the following lemma. 

Lemma 3.2. There exists a unique solution Vp to the problem (j3.7p - (j3.8p on [0,T] x M± x 
M_|_, for any T > 0. Furthermore, for any p > 1, there is some positive constant Cq that 
depends on p and T such that 

||5^yp||ioo(o,T;LP(R±;Tyl.P(R+))) < Co\\d^Voi- ,0)\\ LP fc = 0,l,2, (3.9) 

and the jump of discontinuity [Vp]|^_g satisfies 

\\[Vp]\^^o\\l^(o,T;W^-p(r+)) < Co\[Mx,0)]u=o\- (3-10) 
Proof. Similarly as done for Up, the integral representation for Vp easily yields 

||a^yp(x)||ioo(o,T;H/i.p(R+)) < Co\d^vo{x,0)l k = 0,1,2, (3.11) 

for each nonzero x E M. Taking the U' norm of this inequality in x gives p.9p at once. 
The estimate for the jump of discontinuity of Vp follows similarly by noting that the jump 
[^]|^=o satisfies the similar integral representation to that of Vp. 

□ 

3.1.3 Construction of Vkh 

Similarly, plugging the Ansatz (13. ip into the system ()2.ip yields the profile equation for 
Vkh- 

dtVKH = il + \d,iP\^)dlVKH, [Vkh]\^^o = -i^o]x=o, [dxVKH]u^,=0. (3.12) 
We choose the initial data: 

VKH\t=o{X,z)=T^^^^e^'', ±X>0. (3.13) 

We will derive necessary estimates for the profile Vkh- It turns out convenient to introduce 
a change of variables: 

X = X, i= [\l + \dMs,z)\^)ds, 
Jo 

and write 

VKH{t, X, z) = Vkh {i{t, X,z),X {t, X, z),z). 
In these new variables, we then have 

diVKH = d\VKH. Wh\\^^^ = -Nx=o, [dxVKH]l-^^ = 0, (3.14) 
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with initial data 

VKH\i=, = ^l!^e^^, ±X > 0. (3.15) 

The systems ()3.14p and ()3.15p are the heat equations on each half lines X < and X > 0, 
with z being a parameter. Thus, the Duhamel principle for the heat equation yields a 
candidate for VKH{i,X,z) as 

^ Jo Jo 

for X > 0, and 

VKH{i,X,z)J^e^+^^ f f g{i-s,X;X')e^'dX'ds, 
with the Green function for the heat equation defined by 

Q{1, X- X') = G{t, X-X')- G{t, X + X'), G{t, X) = -^e'^'/^K 

v47rt 

It is straightforward to check that these definitions of Vkh on M.± x M+ indeed satisfy 
the boundary and jump conditions from (|3.14p and (|3.15p . 

Furthermore, similarly to those estimates obtained for Up and Vp, we can easily obtain 

||92Vft-H(2;)|licx,(o,T;iyi.p(R±)) ^ Co\[d^voix,z)]i^^J, ,/c = 0,l,2, 

for each z S M_|_ and for some positive constant Co that depends only on p and T. Going 
back to the original coordinates {t,X), we have thus shown 

||9^VK//(^)||L-(o,T;iyi.f(R±)) < Co\[d%{x,z)]\^^J, (3.16) 
for k = 0,1,2 and for each z G M+. 

Collecting these information, we obtain the following lemma. 

Lemma 3.3. There exists a unique solution Vkh to the problem p.l2p and (j3.13p on 
[0,T] X M-t X for any T > 0. Furthermore, for any p > 1, there is some positive 
constant Cq that depends on p and T such that 

\\dz'^KH\\L^(0,T;LP{R+;W^,P{R±))) < C'o || [S^t'O (^J, •)] U=0 II (K+) ' (3-17) 

fork = 0,1,2. 

Proof. The estimate (|3.17p is obtained easily by taking to both sides of (j3.16p the usual 
norm in z and using the triangle inequality. □ 



10 



3.1.4 Construction of Vb 

Finally, in the "box" where the interactions take place, we obtain from the equation for 
with the Ansatz (j3.ip the following equation for the interaction profile Vb{t,X, Z): 

dtVb = A%Vb, A% = 5i + {dz - ^,^P\^^^^x?, (3.18) 

where 

with the boundary and jump conditions: 

H|,.o = -yKHU^o^ [H]|,.o = -[Vp]u^o^ [5xH]|^.o = 0> (3-19) 
and Vb — )• as X — )• ±oo or Z — )• +oo. Next, we choose the initial data for VJ,: 

which satisfy all the conditions in (|3.19p . thanks to (I3.8P and to (|3.13p . 

We observe at once that these boundary and jump conditions in (j3.2p . (j3.7p . (|3.12p . 
and (13.191) make the Ansatz v^pp defined as in (13. ip smooths out the inviscid solution v*^ 
(at least with regularity) as well as satisfy the correct no-slip boundary conditions (|1.6p 
for the depleted Navier-Stokes system. We will see in the next section that these Ansatz 
v^pp indeed provide a good approximation for v"^, and are sufficient to show the desired 
convergence. 

We will show in this section that the profile Vb exists and we then derive necessary 
estimates to carry out the convergence stated in the main Theorem II. 1[ In fact, we could 
continue our study by employing the Green function of the heat equation on the half-space 
as done previously on the half-line. However, we choose to proceed the analysis by energy 
estimates, as it appears natural for the proof of our desired convergence. 



To begin, it appears convenient to introduce w through 

Vb = 

The function w then solves 



w + l[Vp]\^^^e-\''\, X<0. 



(3.21) 



dtw = AfP^w + J+e-l^l , X >0, 

with boundary and jump conditions: 

[w]\x=o = 0' [dxw]\^^^ = 0, lim w = lim w = 0, (3.22) 

|A|— >-oo Z— >oo 
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and 



w 



|z=o 



Moreover w vanishes at the initial time: 



(3.23) 



w\t=o = 0. (3.24) 

Here, J± in (|3.21|) collects the terms involving the jumps of discontinuity. Direct calculation 
together with a use of ()3.7p gives 



(3.25) 



■^((1 + \d.i^uJ')[VpU, + 2a,V'|^^jazFp]|_o), 
J- = ^((i + |5.V'u.ol')[^p]u^o-29.V'u.o[5^^p]u^o)- 

By applying the estimate (j3.10p obtained in Lemma 13.21 we then have 

\\J±mi,.,v dt < Co|N(x,0)]|_J, (3.26) 

for some Cq that depends on p and T. 

We are able to provide the following estimates. 

Lemma 3.4. There exists a unique solution w to the problem (I3.2ip - ()3.23p on [0, T] x M± x 
M+, for any T > 0. Furthermore, for any p > 1, there is some positive constant Cq that 
depends on p and T such that there holds 

^\ ip+ / |^Z;|^'-2|V^°2^D|2 / \dzwY'-^\^'^j^ zdzw\'^ dXdZ 

< Co(l + ||J±||J^i,. + ||«)||^p i„ 



d „ , 
-r hi' 



Here, V^"^ := {dx,dz - dz'4'\^^odx)- 

Proof. By multiplying by |i()|''~^zZ; to the equation (j3.2ip and integrating it over M^, one 
has 

For the term involving J±, the standard Holder's inequality gives 
I e"^^ J±\wY'-^ dXdZ < \\wf,^^ \\J±\ 
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Here and in what follows, by / < 5 we always mean that / < Cqq, for some positive constant 
Co that only depends on p and T. Now, integration by parts yields 



^ A%w\w\P-^w dXdZ 



(p-l) [ \w\P-^\V%w\^ dXdZ - [{dz- d,ij\^_^dx)w\wr^w\^_^ dX 



dZ 

in which the last term on the right-hand side vanishes due to the jump conditions ()3.22p on 
w and on dxw. Thus, we obtain 

/ A^^^iirZ; dXdZ = - {p - I) [ \w\P-^\V%w\'^ dXdZ - [ dzw\w\^~'^w\^^^ dX. 



(3.28) 



Collecting, we have shown 

+P{P-^) ! \wr''\vfzw\'' dXdZ 

^ ll^ll^p^ II«^±|Il^^ - / dzw\w\^~'^w\ dX. 
For the boundary term, the Young's inequality yields 

\dzw\wf~'^wi^^o\ dX < \\w\^^XlI + \\9zwi^=o\\1^^- 

JR ^ ^ 

The first boundary term can be easily treated by the trace inequality. We treat the second 
boundary term by the energy estimate. To this end, we take Z-derivative of the equation 
(|3.2ip and multiply by \dzw\P~'^dzw to the resulting equation. We simply get 

-^Jdzwfr. = [ (A%dzw + dzJ±e^'')\dzwr^dzwdXdZ. 
Again, by applying the Holder's inequality to the last term on the right-hand side, we have 



^xz 



[ dzJ±e^^\dzw\'"^dzwdXdZ < \\dzJ±\\Lr \\dz 

Next, the integration by parts yields 

/ A%\dzw\^'''^dzwdzw dXdZ 
JrI 

= -{p-l) I \dzwr^\V%dzw\^ dXdZ - [{dz- d,^l^\^^^dx)dzw\dzw\P-^dzw\,^, dX 

JR\ JR 

- j {[d\zMdzM^~'^dzw - d:,i)\^^^{dz - dzi)\,^odx)dz^dz^^^^ dZ 
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in which again the last term on the right-hand side vanishes due to the jump condition 
(j3.22p . By using the equation for w, we can write the boundary term as 

JR 

= - I {{dz - dz^\^^^dxfw + dzil)\,^^{dz - dzil)\,^^,dx)dxyj) \dzw\''~^dzw\z^o dX 
= - [ (dtw - (1 + \d,i^\^_f)dj,w - J±eT^ + d,i^\^_^dj,zw)\dzwr^dzw\^_, dX, 

in which the integral term involving d'j^^w vanishes due to the jump condition p.22p and 
the fact that it is a perfect derivative in X. For the other terms, we note that at Z = we 
have 

dtw - (1 + \d,iP\^J^)dj,w - J±e^^ = dtVKH - (1 + \d,iP\^J^)dj,VKH + c±e=F^ 

for some constant c±; here, the last identity was due to a use of the equation for Vkh- 
Thus, using this and the Sobolev embedding, we have 

- (dz - dzTp\^^,,dx)dzw\9zw\^~'^dzw\z^o dX 
Jr 

< \\dzw\^^o\\LPj, 

1/p 



< 



\dzw\\L^^^ + ( [ laz^r^Vl^yp dZdX 



Thus, applying the Young's inequality to the last term and combing all the above estimates, 
we obtain 

^\\dzw\\lv + I IdzwlP-^lV^j^'z^zw]^ dXdZ 

dt XZ Jj^2^ 

< 1 + WOzwWlp^ + \\dzJ±\\ LP \\dzw f]~p^ ■ 
This together with the estimate ()3.28p yields the lemma at once. □ 

We also obtain the following X-derivative estimates. 

Lemma 3.5. For any solutions w to ()3.2ip - ()3.23p . there holds 

i\\dxw\\%^r..+ f \dxwr^\vt'zdxw\'dXdZ+ [ {d'xzwr'l'^'^x z^xz^l' dXdZ 

< i + iiJ±r^,. + ii«)r^,.. 
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Proof. The proof of this lemma fohows word by word of that of the Lemma 13.41 upon 
noting that the jump of discontinuity of dj^w across {X = 0} can be computed through the 
equation (j3.2ip for w to give = — J^. Note also that we may need to apply 

the Sobolev embedding: 



i/p 



We thus omit the further detail of the proof of the lemma. □ 

To conclude this subsection, we summarize our estimate for Vh in the following lemma. 

Lemma 3.6. There exists a unique solution Vh of the problem ()3.18p with the boundary and 
jump conditions (I3.19P and initial data (13.20p . Furthermore, for any p > 1, there exists 
some positive constant Cq that depends on p and T such that there holds 

sup \mtW,,, + r [ dXdZdt 
0<t<T ^xz Jq Jir2 (3.29) 

< Co (l [vo{x, 0)]|^^„ \P + II [vo{x, Olu^o \\w^,p(R+) + ll^o(-, 0) ||^i,p(jj^)) . 

Proof. This is a collection of estimates from Lemmas 13. 41 and 13. 5| the estimate (j3.26p on J±, 
the jump estimate (j3.10p from Lemma 13.21 and a use of the standard Gronwall inequality. 
Indeed, Lemmas 13.41 and 13.51 inparticular yields 

sviP \\w\\Ia P + [ [ \dxw\^~^\dlw\'^ dXdZdt 

0<t<T ^xz Jo Jm?^ 

is bounded. This together with the standard Young's inequality yields that 

/ / \dlw\P dXdZdt = [ [ Idxwl^^ldxwl'^PldlwlP dXdZdt 
Jo Jr^ Jo Jr^ 

\dxw\f dXdZdt+ [ [ \dxw\P-^\djcw\'^ dXdZdt 
I Jo Jml 



< 







is also bounded. The lemma is proved. □ 
3.2 Remainders 

We observe that, with the above profiles, the functions {ul^pp,v^pp) given by the formula 
(j3.ip satisfy 

dtKpp = eS^n^pp + 
dtvlpp + {ulpp - uo)d^vlpp = eAtXapp + ^^ 

«PP><pp)u=o = 0' (3-30) 

[9x<pp]|,^o = 0, 
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where direct computations give 

= edluo, 

and 

+ e(l + |a,V'|'))92yp - ^^ed.^dl^Vp - edl^d^Yp - V^d^v^^^. 

Remark 3.7. We note that E"" contains a singular term: -^^^i^^zj \ft)dxVKH^-,xl \ft^z) 

in Updxvl^pp. This singular term is a-priori not better than bounded in and thus we can't 
obtain the convergence in the space this way. However, its LP norm has an order of 
uuJiich tends to zero in L'P, for 1 < p < 2, as e — )• 0. 



Proposition 3.8. For all p > 1, there hold uniform estimates: 

l|i^"(i)llLi < e||9^o||Li, r WE^mi. dt < an{p,T)e'-P/^ (3.31) 

Jo 

for some positive constant Cin{p,T) that depends continuously on the initial data, the dis- 
continuous jump of the Euler flow vq, the number p > 1, and the time T. More precisely, 
the constant Cin{p,T) is bounded by 

C7o(|N(x,0)]|_JP + ||[7;o(x,-)]|_ollwi.f(R+) + ll^o(-,0)||^i,,(jg^) + l|At;o||^p^(K^^K+)) 

for some Cq that depends only on T and p. 

We now give a proof of the Proposition [3]8l The first estimate is clear from the definition 
E^ = ed'^UQ. We prove the second estimate. We will use the following simple lemma. 

Lemma 3.9. For any reasonable function u = u{z), there holds 

\\u{z/^^e)\\Ll{w.+) < e^^'^lk(^)llL|(K+)- 

Proof. It is clear by changing of variable from z to z/y/e. □ 

We now check term by term in E^ . The term eA^^uo is clear, giving the contribution 
of e||A^;o||^P^(R^xR_^). Next, note that 

Thus, the estimate ()3.29p for Vi, precisely gives us the desired estimate for (A^^ — 
A^°2)Vfe, after a change of variables (x, z) to {X, Z) with X = x/ ^/e, Z = z/i/e, yielding a 
small factor of e. Similarly, for all the terms: 

-2V^a,V'92^VKiT-^/^a,V5xV'xH+e9,V;^H+e(l+|5,V'|2))a2yp-2V^a,V'a2^yp-ea,V5x^P 
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the estimates from Lemmas 13.11 13.21 and 13.31 immediately yield that the norm of these 
are bounded by 



Finally, let us treat the term Updxvl^pp. From the definition of v^pp, the singular terms 
in Updxvlpp are 



^Up{t, Z)dxVKHit, X, z) + ^Up{t, Z)dxVb{t, X, Z). 
Ve ye 



We then use the Lemma 13.91 to treat these singular terms. For example, we compute 



/ \Up{t,Z)dxVb{t,X,Z)\P dxdz < /" \Up{t,Z)dxVbit,X,Z)\P dXdZ 

< e'-^'^'lKWUldxV.Wl, 

z ^XZ 

< ei-p/2 

Other terms are entirely similar. This completes the proof of the estimate (j3.3ip . and thus 
the Proposition! 



3.3 Convergence 

We are ready to prove the convergence stated in the main theorem. 
Now we consider the solutions i?", K" of the following problem: 

dtR" + {Up + W')dxW = eAtR" - R'^d.vlpp + 
(i?",i?-)|^^„=0, 
lim =0, 

z— >+oo 



(3.32) 



Then the functions {u^,v'') defined by 



u'it,z) = ulpp{t,z) + R''{t,z) 
v^{t, X, z) = vlpp{t, X, z) + R''{t, X, z), 



satisfy the equations (|2l41 -(l2:2])-(l2:ill. 

The well-posedness of the problem (j3.32p follows at once from the following a-priori 
estimates: 
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Lemma 3.10. There hold 

d 
di 



,11 IIP— 1 



I'l^'llk. \RT-'\^tli"\' dxdz < (||i?"«ppllLL + 

where vt,z ■= {dx,dz - dzi>dx)- 

Proof. Multiply by ii" and \R^Y'~'^R^ the respective equations in (j3.32|) and integrate the 
resulting equations over or R^. The claimed estimate for i?" is straightforward. For 
the estimate, we have 

-4/ \R''F+ f (Ur, + R'')\R''F~^dxR''R^ dxdz 



p dt 



= J ^ [eAtR" - R^'dxvlpp - E''^ liJ^I^-^i?'' dxdz. 

We first note that integration by parts yields 

/ {Up + R''')\R''\P-'^dxR^R''' dxdz = [ {U„ + R'')dx(^-^-^) dxdz = 

jRl Jul \ p J 



and 



:/ AtR''\R''\P-^R'' dxdz = -e [ iR^l^-^lVt .R"]'^ dxdz, 



upon noting that there is no contribution on the boundary and the interface due to the 
vanishing boundary and jump conditions for R". Finally, the standard Holder inequality 
yields 

/ {\R-d,vl^^\ + \E^\)\R^rUxdz < (\\R-d,vlp^h.^ + \\E^hAmi^'- 
Collecting these estimates together proves the lemma. □ 

It is straightforward to verify that 

{ulpp{t,z),vlpp(t,x,z)) ^ v" in L°°(0,T;L2(M+) x LP{R x M+)), 

as e — 0, upon using the estimates on the profiles from Lemmas 13. H 13. 2|, 13. 3|, and 13.61 and 
the Lemma 13.91 Therefore in order to prove Theorem I l.H it remains to prove that 

{R''{t,z),R''{t,x,z)) ^0 in L°°(0,T;L2(M+) X X (3.33) 

as e — )• 0. 
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From Lemma [3 . 1 1 and Proposition l3.8l by the standard ODE estimate and the Gronwall 
inequahty, we immediately obtain uniform bounds 

with noting that ||-E"||^2^2 < e and = 0. In addition, this estimate yields 



Jo * " 



u\\2 < J 



This together with the bound H^i^^^appllioo^^p ^ which again follows from the estimates 
on the profiles, yields 



In addition, the second estimate from Lemma 13.101 implies 



which gives 



xz xz V 



The Gronwall inequality then yields 



\\J}V\\P < 3/4 , l-p/2 

11-" IIl-lL ^ e +e 

which tends to zero as e — )• 0, for p <2. 

This ends the proof of the convergence (j3.33p , and thus of Theorem II. li 
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